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We investigate the Kramers- Wannier approximation for the three-dimensional (3D) Ising 
model. The variational state is represented by an effective 2D Ising model, which contains 
two variational parameters. We numerically calculate the variational partition function using 
the corner transfer matrix renormalization group (CTMRG) method, and find its maximum 
with respect to the variational parameters. The calculated transition point K c — 0.2184 
is only 1.5% less than the true K c ; the result is better than that obtained by the corner 
transfer tensor renormalization group (CTTRG) approach. The calculated phase transition 
is mean-field like. 

§1. Introduction 

In 1941 Kramers and Wannier ^ proposed a variational approximation for the 
two-dimensional (2D) Ising model, which is called 'Kramers- Wannier (KW) approx- 
imation' today. The feature of the approximation is that the variational state is 
constructed as the thermal equilibrium state of the ID Ising model in an effective 
magnetic field. From the modern view point, their variational state can be regarded 
as an example of so called the matrix product state. 3 )' 4 ) More than 20 years later, 
Baxter improved the KW approximation by introducing additional degrees of free- 
dom into the variational state; he reformulated the variational principle of the KW 
approximation using the corner transfer matrix (CTM). 5 )~ 7 ) It has been known that 
the variational property in both the KW approximation and Baxter's CTM formu- 
lation has many aspects in common with that in the density matrix renormalization 
group (DMRG) 8 ) " n ) and the recurrent variational ansatz. n )' 12 ) 

The transition temperature and the specific heat of the 2D Ising model calculated 
by the KW approximation are more accurate than those obtained by the mean- 
field approximation and the Bethe approximation. 2 ) It is expected that the KW 
approximation is a good non-perturbative method also in higher dimensions. We 
therefore investigate the KW approximation for the 3D Ising model by way of the 
maximization of the Rayleigh ratio 

A _ g) rin 

WT' ( } 

where T is the 'layer-to-layer' transfer matrix, and \V) is the variational state rep- 
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resented as a 2D generalization of the matrix product state. Though the varia- 
tional formulation is quite simple, such a generalization of the KW approximation 
to 3D systems has not been investigated so far. This is partially because there is 
no analytical tool to calculate (V|T|V) and (V|V) in eq. (IT) that are partition 
functions of unsolvable 2D lattice models. In this paper we overcome the problem 
by numerically calculating A using the corner transfer matrix renormalization group 
(CTMRG), 16 )> 17 ) which is a variant of the DMRG for 2D classical systems n )' 18 ) " 20 ) 
formulated via Baxter's CTM. 5 ) " 7 ) Our approach shown in the following can be re- 
garded as the DMRG applied to 3D classical systems. It should be noted that the 
formulation of the KW approximation in eq. (IT) is related to the tensor product 
variational formulation for 2D quantum systems. 13 )~ 15 ) 

In the next section, we introduce the KW variational state \V) for the 3D Ising 
model, and present the concrete definition of eq. (IT). In §3, we explain the way to 
apply CTMRG to the variational formulation, and then show the calculated spon- 
taneous magnetization and the internal energy. Conclusions are summarized in §4, 
and we discuss several possible improvements to the formulation of the KW approx- 
imation in 3D. 

§2. Variational Formulation in 3D 

We consider the 3D Ising model on a simple cubic lattice of the size N x N x L 
in the X, Y, and Z direction, respectively, where on each lattice point (i,j,k) - 
the position in the fc-th spin layer — there is an Ising spin afj = ±1. We 

assume open boundary conditions in both X and Y directions, and periodic boundary 
condition in Z-direction. The Hamiltonian of the 3D Ising model is 

H = -JB^j + offotj + ogo*-) , (2T) 

ijk 

where we have used the notation %' = i + 1, j' = j + 1, and k' = k + 1 for book 
keeping. Throughout this paper we consider the ferromagnetic case J > 0. The 
partition function of the system is expressed as 

Z = ]Texp(-/3tf) = TrT L , (2-2) 
M 

where T is the layer-to-layer transfer matrix, and the sum is taken over for all the 
spin configurations. In the following, we consider the symmetrized transfer matrix 
T(= T T ) constructed as a product of local Boltzmann weights 

T(a k '\a k )=l[W l k j , (2-3) 

ij 

where a k and a k ' represent spin configurations in fc-th and fc+l-th layer, respec- 
tively, 21 ) and W k j is the local Boltzmann weight defined by 

W k j = ex P j^j" ( a i'j a ij + a ij' a ij + a i'j' a i'j + a i>j' a ij> 



Kramers- Wannier Approximation for the 3D Ising Model 



3 



i j i j a i'j a i'j °i j /(J i j' + a i'j' a i>j' 
+ 0-%0-\j + Crf j'ffj + O-i'fCTi'j + <Ti'j'(Tij') | , (2-4) 

with K = (3J. In the r.h.s of the eq. (2-4), the 12 terms correspond to the 12 edges 
of a local cube, and the coefficient 1/4 of K denotes that each bond between the 
nearest neighbor spins is shared by 4 adjacent cubes. 

As an introduction to the KW approximation for the 3D Ising model, let us 
consider a special mean-field approximation, which replaces all the Ising spins of - 
except at the fc-th spin layer a k by their expectation value (a). The approximation 
draws the effective Hamiltonian for the k-th spin layer 

H(a k ) = -jJ2(4 3 o-i 3 + oiyo-^j + 2{a)a k ij ) , (2-5) 

which is nothing but the Hamiltonian of the 2D Ising model under the mean field 
2J(a) imposed from both the up and down sides of the fc-th layer. In this mean- 
field framework, the spin profile in the k-th layer is given by the weight P(a k ) = 
exp{— (3H(a k )}. It is expected that the mean field weight P(a k ) well approximates 
the appearance probability of the layer-spin configuration, and that its square root 

y/P(<r k ) = II ex Pi ^X^fj + + °ij> + Oi>j>) 

ij K 

K k k \ 

+ -ji^i'j^ij + a ij' a ij + a i'j' a i'j + a i'j' a ij') f ) ( 2 ' 6 ) 

which is proportional to T(a k \ (a)), can be used for the variational state in eq. (IT). 

Such a direct usage of \J P{o k ) as the variational state, however, has a short- 
coming in the paramagnetic region, where {a) is zero and yj P{a k ) has no adjustable 
parameter. Following Kramers and Wannier, we introduce an additional parameter 
to the nearest neighbor coupling term in eq. (2-6). The variational state (in the 
product form) is then given by 

V(a k ) = l[U k j , (2-7) 

ij 

where the local factor U k , is defined as 

1 j 




with two variational parameters h (= effective magnetic field) and g (= effective 
nearest neighbor coupling). The variational state V(a k ) in eq. (2-7) has at least one 
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variational parameter g even when the system is paramagnetic (h = 0). Substituting 
T(a k | a k ) and V(a k ) to the variational formulation in eq. (1-1), the Rayleigh ratio 
- the approximation partition function per layer of the size N x N — is expressed 
as 

E {ak ' } Wk} V(v k ')T(a k '\a k ) V(a k ) 
where the denominator of the r.h.s. 



= - "v . (^) 



{cr fc } ij {cr fe } ij I 

+ 7;( a i'j a ij + a ij ,(7 ij + °i>j>°i'j + 

is a partition function of an effective 2D Ising model parameterized by h and g. 
Similarly, the numerator 

B N = V(a k ')T(a k '\a k )V(a k )= £ U^Wi (2-11) 

{<T fc '}{<7 fe } {^'}{^} V 

is a partition function of a two-layer Ising model parameterized by h, g, and K. 

§3. Numerical Result 

The goal of the KW approximation is to find out the pair of h and g — as 
functions of K — that maximizes the variational partition function per site in the 
thermodynamic limit: 

z(K,h,g) = lim (Xn) 1 ^ 2 = lim (Bn/An) 1 '"* . (3-1) 

TV— >oo N—*oo 

In order to find out the maximum of z(K,h,g) in the h-g parameter space, we 
calculate z(K,h,g) for various values of h and g via the numerical calculation of 
An and for N = 3, 5, 7, . . . up to a sufficiently large N. After that we search 
the maximum of z(K, h,g) in the h-g plane. From the numerical point of view, it is 
better to use the formulation 

z(K h g) = lim ( Bn +^ n + 2Bn \ ^ (3.2) 

in order to accelerate the numerical convergence with respect to TV, rather than just 
taking the limit TV — ► 00 directly to (Bn/An) 1 ^ . 

We use the CTMRG method 16 )> 17 ) for the calculation of and Bn, since the 
method enables us to obtain An and B jy very rapidly and accurately. We keep 
m = 32 states in the CTMRG calculations, and obtain An and B n up to N = 150; 
the condition is sufficient for the precise determination of z(K,g,h). Throughout 
this section we set J = 1 and thus K = [3. 
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Fig. 1. The internal energy £ as a function 
of K. The transition point is observed at 
K c = 0.2184. 



M 



1.0 
0.8 
0.6 

1 

0.4 
0.2 



0.0 



0.2 



0.3 

K 



0.4 



Fig. 2. The magnetization M as a function of 
K. The solid line denotes eq.(3-4) obtained 
by Monte Carlo simulations. 22 '' 23 - 1 In the 
very vicinity of the calculated transition 
point, M is proportional to \JK — K c . 



First, let us see the calculated results of thermodynamic quantities. Figure 1 
shows the internal energy per cube, which is given by 

E = -{°ij°v j) ~ (o'ijO'ij') ~ (vijVij) , (3-3) 

according to the formal thermodynamic relation E = — In z(K,h, g). We have 
checked that the numerically calculated In z(K, h,g) really satisfies the relation. The 
plotted data in Fig. 1 has a kink at K c = 0.2184, which is the transition point from 
the paramagnetic state to the ferromagnetic state. The calculated K c is about 
1.5% smaller than one of the reliable critical point K^ lc = 0.2216544 ± 0.000005 
determined by Monte Carlo simulations. 22 )' 23 ) The discrepancy 1.5% shows that the 
KW approximation for the 3D Ising model is more accurate than that for the 2D 
Ising model; for the latter, the critical point K c = 0.4122 calculated by the KW 
approximation is 6.5% smaller than the exact one K c = 0.4407. 24 ) 

In Fig. 2, we draw the spontaneous magnetization M = (o~fj)- For comparison, 
we also show Tarpov and Blote's Monte Carlo result for the cubic lattice containing 
up to 256 3 spins, (see eq. (10) in Ref.[22]) 

m mc = t o.32694i09( L6919045 _ o.34357731t a50842026 - 04 2 5 7 2 3 66t) , (3-4) 

where t = 1 — K^ lc /K, and this expresson is very accurate in t < 0.26. The KW 
results agrees with eq. (3-4) in almost whole the region of K. In the very vicinity 
of the transition point K c < K < K c + 0.01, the calculated magnetization deviates 
from eq. (3-4) and behaves as 5.171y/K — K c approximately. In principle, the phase 
transition observed by the KW approximation is mean-field like. 25 ^ 

Let us see the properties of the variational parameters g and h, respectively, 
since they are closely related to E and M. As the internal energy E in Fig. 1, the 
parameter g shown in Fig. 3 has a kink at the calculated K c , and is always larger than 
K. The parameter h shown in Fig. 4 is approximately 0. 5516a/ K — K c in the vicinity 
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Fig. 3. The variational parameter g as a func- 
tion of K. 
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Fig. 4. The variational parameter h as a func- 
tion of K. In the vicinity of the transition 
point, h is proportional to \/ K — K c . 



of the calculated i^ c ; 2/i ~ KM is approximately satisfied in the neighborhood of 
K c . 

§4. Conclusion and Discussion 

We have applied the KW approximation to the 3D Ising model, representing 
the variational state as the thermal equilibrium state of an effective 2D Ising model. 
We have calculated the variational partition function numerically using CTMRG 
method, and maximized the function with respect to the variational parameters g 
and h. 

The KW approximation draws the spontaneous magnetization fairly well in wide 
region of temperature, compared with the Monte Carlo simulations. The calculated 
transition point K c = 0.2184 is only 1.5% smaller than one of the most reliable 
K c determined by Monte Carlo simulations; 22 )' 23 ) the K c obtained in the KW ap- 
proximation is better than that obtained by the corner tensor renormalization group 
(CTTRG). 26 ) The critical behavior observed by the KW approximation is mean- field 
like in the very vicinity of the transition point. 

There are at least two ways to improve the variational state used in the KW 
approximation for the 3D Ising model. A way is to introduce additional variational 
parameters into the trial state V(a k ). Thermal equilibrium state of arbitrary 2D 
classical lattice models, such as the multi-layer 2D Ising model and the Ising model 
with next nearest neighbor interactions, can be the candidates for V(a k ). It is 
straightforward to apply the CTMRG to such a variational state to evaluate the 
variational partition function. However, the number of variational parameters is 
limited by the numerical effort to find out optimal variational parameter sets. The 
other way is to introduce block spin variables into each local factors in V(o k ). 15 -*' 19 -* 
Although this approach contains much more variational parameters than the former 
way, we can treat the optimization problem more systematically as was done in the 
CTTRG. 26 ) In both of these improvements, the key point is to find out the best 
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variational parameter set quickly. 

We finally comment that the formulation of the KW approximation for the 3D 
Ising model presented here can be applied to various 2D quantum spin systems. The 
generalization is simply to replace the 2D tensor product state in Hieida's DMRG 
formulation 14 ) by the thermal equilibrium state of 2D classical lattice models. This 
is not a trivial simplification, since the relation between the 2D KW variational state 
and the 2D tensor product state have not been clarified yet, unlike the trivial relation 
in ID. 
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